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ABSTRACT 


Laboratory tests conducted to determine the effects of wave groups on rubble 
mound breakwater stability are analyzed. A literature review outlining the controversy 
regarding whether or not wave groups effect breakwater stability 1s given. Resolution of 
this controversy is proposed using two independent parameters that characterize wave 
groups. The two parameters are an envelope exceedance coefficient œ and spectral 
shape y. The incident time series must be known to compute a, and an algorithm is 
given that will resolve incident and reflected wave time series from closely spaced wave 
gauges. The digital to analog simulation of laboratory waves with specified wave group 
characteristics is reviewed. It is recommended that further studies be conducted to 


parameterize the envelope exceedance coefficient a with a non-biased groupiness factor. 
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INTRODUCTION 


The design of rubble mound breakwaters is one of the oldest and most common 
applications of coastal engineering methodologies. Design of rubble mound structures is 
most commonly based on the empirical methods presented in the Shore Protection 
Manual (SPM 1984). Due to some failures of structures designed with the SPM 
methodology, an effort to improve the SPM design methodology is undertaken. The 
design methodology presented in the SPM is based on empirical formulas developed 
from experiments with monochromatic waves. Most of the more recently proposed 
design criteria are based on laboratory experiments using irregular waves. The effects of 
spectral shape and wave grouping will be targeted. As a proposed improvement, factors 
such as permeability, wave period, and storm duration may have secondary effects on 
rubble mound breakwater stability. 

It has been well-documented that waves in random seas tend to form groups not 
accounted for by the purely Gaussian random wave model. Many of the coastal 
engineering problems associated with wave groups have been documented by Medina 
and Hudspeth (1990). In the last two decades, controversy has arisen as to whether 
spectra composed of grouped waves cause different levels of damage to rubble mound 
Structures than spectra of ungrouped waves. Medina et al. (1990) have found that most 
of the controversy may be resolved if two independent spectral related parameters are 
chosen to characterize the wave groups (Mase and Iwagaki 1986). One parameter 


characterizes the length of wave groups and the other characterizes the magnitude of the 
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energy flux in a group of waves. 

Intermediate scale laboratory studies to determine if wave groups affect rubble 
mound breakwater stability were performed in the 2-dimensional wave channel at the O. 
H. Hindsdale Wave Research Facility (OHH-WRF) at Oregon State University. In 
order to determine the effects of wave groups on rubble mound structures, sea states 
with identica! amplitude spectra but different wave group characteristics were simulated. 
The resolution of incident and reflected wave time series from the random wave 
simulation is required in order to evaluate the effects of wave groups. In order to 
quantify the two wave group characteristics proposed by Mase and Iwagaki (1986), the 
incident time series must be known. Medina et al. (1994) used a method developed by 
Kimura (1985) and Fassardi (1993) that extends an algorithm proposed by Goda and 
Suzuki (1976) to resolve the incident and reflected wave time series from closely spaced 


wave gauges. 





CHAPTER I: Correlation of Wave Groups with Armor Damage 


Carstens, et al. (1967), first observed the effect of wave groups on breakwater 
damage when he analyzed two time series with significantly different groupiness, but 
similar spectra and wave statistics. The time series with the more grouped waves, taken 
from a field site 1n the Barents Sea, created. more damage than the time series generated 
from the theoretical Neumann-spectrum. Johnson, et al. (1978), as a result of their 
laboratory studies on the effects of wave grouping on breakwater stability, concluded that 
wave grouping was an essential parameter in the model testing of rubble mound 
breakwater stability. Wave trains with similar spectra and wave statistics, but different 
eroupiness, caused significantly different levels of damage, the greater damage being 
caused by the grouped wave trains. Burcharth (1979) also concluded that wave grouping, 
or the succession of the waves, 15 Important ın model design. 

Since the energy spectrum does not contain the phases of the individual wave 
components in the sea surface 1t is necessary to also resolve the phases of the individual 
wave components to reproduce an accurate succession of wave heights (Rye 1982). In 
his summary of breakwater design, Bruun (1989) recorded that the resolution of the 
phase spectrum identified apparent order where none is assumed in the Gaussian 
random wave model. He concluded that these ordered wave groups are more damaging 
than other waves. 

In contrast, van der Meer (1988) found that the spectral shape and grouping of 


waves had little or no influence on the stability of breakwater armor layers provided that 





the average spectral period was used to simulate the sea state rather than the peak 
period. Hall (1994), in a study on bermed breakwaters, concluded that wave groups are 
only important during the stage of incipient motion, and that once motion 1s initiated 
wave groupiness is no longer important. 

Medina et al. (1990) found that much of the controversy could be resolved if the 
two independent parameters of wave groups proposed by Mase and Iwagaki (1986) were 
considered. One parameter measures the run length or the number of high waves in a 
time series, and a second parameter quantifies the magnitude of the variation of wave 
energy. In their laboratory tests, Medina et al. (1990) tested four different wave height 
time series proposed by Mase and Iwagaki (1986) shown in Figure 1.1, while keeping 


other secondary factors that affect breakwater stability constant. 
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Figure 1.1 - Classification of wave height time series (Fassardi 1993) 





Figure 1.1 shows that identical spectra may have different magnitudes of wave height 
variability, which may cause different levels of damage. Likewise, different spectra may 
have the same magnitude of wave height variability, and thus cause the same level of 
damage. Medina et al. (1990) simulated truncated Goda-JONSWAP spectra with 
different run lengths and different magnitudes of variation of wave energy according to 
Appendix B. 

Medina et al. (1990) concluded that the spectral shape and run length alone were 
not sufficient parameters to accurately describe wave groupiness, but that another 
parameter that measured the variation of wave energy more completely described the 
wave group. Medina et al. (1990) correlated an envelope exceedance coefficient o that 
measures the variation of wave energy with stability, and demonstrated that random 
wave trains with the same spectral shape may produce different levels of damage. 
Medina et al. (1994) found no significant correlation between the spectral peakedness y 
and armor damage, but concluded that the envelope exceedance coefficient a was 
necessary to evaluate the effects of random waves on structures in both physical and 


prototypical models. 





CHAPTER II: Derivation of Envelope Exceedance Coefficient a 


From their analysis of storm waves, Mase and Iwagaki (1986) recommended using 
at least two parameters to quantify wave groups, one parameter that measures the 
magnitude of a sequence of high waves in a time Series, and a second parameter that 
measures the magnitude of the variation of wave energy. Both parameters are shown 


schematically in Figure 2.1, where A(x,t) is the wave envelope function from the Hilbert 


transform defined by Medina and Hudspeth (1990). 


#ዖ run length is a meaure of the number of high waves in a group 
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Figure 2.1 - Wave group parameters 





The peak enhancement factor y of the Goda-JONSWAP spectrum (Goda 1985) 
controls the sharpness of the spectral peak and was selected to characterize the 
magnitude of the sequence of high waves in the time series. Rye (1980) related the 
magnitude of the sequence of high waves, run length, to the peak enhancement factor y 
of the JONSWAP spectrum. The coefficient a was selected to characterize the 
magnitude of variation of wave energy. The envelope exceedance coefficient o is a 
dimensionless measure of the wave energy exceeding that associated with a prescribed 
DSXUht e. g., Hijo. 

In a two dimensional wave channel, the wave height function H(x,t) is defined by 

H(x,t)= 2A(x,t) (2.1) 
If At = a constant sampling time interval, the discrete wave height function at a fixed 
location may be expressed as H(nAt). Defining a normalized measure of the variation of 


wave height as 


| AH, dd (ከል » 


22 
| H ۱ ) 


1/10 


and the Heaviside step function (Farlow 1982) of AH, as 


0 AH <0 


then the magnitude of the variation of wave energy above H;,,¡p of the time series is given 


as 





መድን 


N 
ois 2 Y (AH, PU(AH,) (2.4) 
N በ=1 


where N is the number of discrete data points in the time series. Normalizing Eq. (2.4) 


by the expected value of a’ yields the envelope exceedance coefficient a defined as 


f 





)2.5( د کن 
] 16 
where E[a’] is defined as‏ 
co H 2‏ 
E[c'] - | - 1 p(H) dH (2.6)‏ 
f DO‏ 


110 


where H is the wave height function at a fixed location H(t), and p(H) 1s the Rayleigh 


probability density function of H(t) (Appendix A) defined by Chakrabarti (1987) as 


pup 








er? 
g ex | 2 | (2.7) 
dm, 8m, 


and m, is the variance of the time series. If the following substitution is made (Sarpkaya 


and Isaacson 1981), 


1 mj ug) 


Eq. (2.6) may be solved numerically to obtain E[o'] —0.001346. 





CHAPTER III: Resolution of Incident and Reflected Wave Time Series 


In order to determine experimental wave heights and to calculate envelope 
exceedance coefficients, the incident time series must be known. Because of reflection 
from models in a 2-dimensional wave channel, a method to resolve the incident and 
reflected random wave time series from wave gauges must be used. Bad data points in 
wave records caused by wave gauge malfunction are smoothed using the FORTRAN 
program CLEANREC that employs a smoothing spline package developed by H. J. 
Woltring (1986). A Fast Fourier Transform (FFT) analysis is used to compute the raw 
wave spectrum from which the incident and reflected random wave time series may be 
determined. Goda and Suzuki (1976) modified a method developed by Thornton and 
Calhoun (1972) to resolve the incident and reflected wave spectra from wave gauges 


spatially separated bya distance Af (see Fig. 3.1). 
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Figure 3.1- Definition Sketch for two wave gauges 





Kimura (1985) and Fassardi (1993) extended the Goda and Suzuki algorithm to 
include wave phases so that incident and reflected time series could be resolved. 

Goda and Suzuki (1976) analyzed two simultaneous wave records recorded by 
closely spaced wave gauges aligned in the direction of wave propagation by an FFT. If 
the incident n and reflected (random wave time series at the i? wave gauge at frequency 
w are given by 

n,= a cos(kx,- wt+ €) (3.1a) 


f= b cos(kx,+ wt+ B) (3. 1b) 


then the wave profile at the i^ wave gauge may be expressed as 


(+ ی(‎ 7 Acoset* Bsinet G2) 
where 
A,= a cosy,+ b cosy, (3.38) 
B= a siny,- b siny} (3.3ከ) 
A,= a cos(kAfl + y )+ b cos(kAf + y) (3.3c) 
B,= a sin(kAf + W)- b sin(kA + y) (3.3d) 


ከግ 


where the spatial phases of the incident and reflected waves, respectively, are defined as 


y, kx,* £ (3.4a) 
j.- lx, B (3.4b) 
x," x,* At (3.4c) 


and x,and x,are the wave gauge positions as shown in Figure 3.1. The coefficients A; 
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and B, may be related to the complex-valued FFT coefficient X by (Dean and Dalrymple 


1984) 


Mee)! i (3.5) 





where j= y- 1 . Goda and Suzuki (1976) solve Eqs. (3.3) simultaneously for estimates 


of the amplitudes a and b listed in Table 3.1. 

For irregular waves, the time series are the superposition of many wave 
components given by Eqs. (3.1). Because Eqs. (3.3) apply to each Fourier component in 
an irregular wave record, linear superposition may be used to resolve time series from a 
spectrum of Fourier components (Goda 1985). Goda and Suzuki (1976) observed that 
spectral estimates diverged (toward infinity) near frequencies satisfying the condition 
5 tor n=0,1,2,..., because the term | sin kA? | in the denominator of the 
equations for the amplitudes a and b in Table 3.1 becomes small and errors from noise 
are amplified. Therefore, the wave gauge spacing Af determined the frequency limits of 
a band pass filter from which the waves could be separated into incident and reflected 
time series. The wave amplitudes a and b could be resolved effectively for the interval 
O.laskAf<s0.97. Goda and Suzuki (1976) recommended the following effective band 


pass limits for experimental conditions: 


005 DE MS, (3.68) 


where L,,, and L,; denote the wavelengths corresponding to the lower (f,,) and upper 


(fmax) frequency limits, respectively, of the band pass filter. Goda and Suzuki (1976) note 
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that the effective range may be taken slightly wider than given in Eq. (3.62) such that 


0.031 € At € 0.45L,.. (3.6b) 


Goda and Suzuki (1976) also recommended that the wave gauges be located at least one 
wavelength away from both the breakwater toe and the wave generator. 

Kimura (1985) extended the Goda and Suzuki algorithm to include situations 
where each incident wave component in the Fourier spectrum results in a different 
reflection coefficient from the structure. Kimura (1985) defined incident and reflected 
wave profiles by Eqs. (3.1), and the spatial phases by Eqs. (3.4a,b). Eqs. (3.1) are similar 
to Kimura’s Eqs. (9 & 12) on pages 62-63 (1985) if the sign of the temporal term is 
changed from minus to plus in his Eq. (12). Kimura (1985) also solved for the spatial 
wave phases defined in Eqs. (3.4a,b). Values for a, b, Y,, and Y may be calculated using 
the equations listed in Table 3.1 for Kimura. To improve the resolution of incident and 


reflected waves, Kimura recommended a slightly more conservative data window of 


QNID. «Af 0.35, (3.7) 


Fassardi (1993) followed the approach of Kimura (1985) to extend the Goda and 
Suzuki algorithm (1976). However, he defined the incident 7 and reflected ¢ time series 


at the i^ wave gauge as 
AMOS (O m= 6) (3.5) 
ከ 5 ር050ሟ(ነ 6! ዐ) (3.8b) 
If the spatial phase of the incident and reflected waves, respectively, are defined as 


y,= kx - € (3.9a) 
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Y= kx, + p (3.9b) 


rather than by Eqs. (3.4a,b), then the amplitudes and spatial phases are given by the 
equations listed in Table 3.1 for Fassardi (1993). The only difference between the 
Fassardi and Kimura algorithm is the sign of €1n the definition of the incident time 
series phase. This sign change results in corresponding sign changes in the Fourier 


coefficients that are now given by 


A,= a cosp,+ b cosy, (3.10a) 
B,= a sind, + تا‎ siny, (3.10) 
A,= a cos(-kAf+ (+ b cos(kAl+ Y.) (3.10c) 
B,= a sin(-kAl+ ሃዕ ፥+ ኩ ቋበርኗሷያ ቶ ህ;ን (3.100) 


These four equations are then solved simultaneously to obtain the equations for a, b, v, 
and vr listed in Table 3.1 for Fassardi (1993). 

In the experiments documented by Medina et al. (1994), the Fassardi algorithm 
was used to resolve the incident and reflected time series. Three sonic wave gauges 
aligned in the direction of wave propagation were centered 10m from the toe of the 
breakwater and spatially separated by A@=1.22m according to Figure 3.2. 

Clean time series from each gauge were used in the FORTRAN program 
RECANSYG to resolve the incident and reflected time series. RECANSYG outputs the 
raw time series and spectra resolved from the wave gauge system in Figure 3.2. The 
program computes the time series mean and variance. As well, it computes the 


reflection coefficient, significant wave height, and mean period from the moments. The 
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Table 3.1- Summary of incident and reflected wave amplitudes and phases 
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incident and reflected time series were resolved using Gauges 1 and 2 and Gauges 2 and 
3 so that the algorithm could be verified by duplication. A sample run is analyzed. As 
explained in APPENDIX B - Simulation of Sea States, Run EIPIL7 represents a sea 
state simulated with the following parameters: 

El = Envelope 1 ; y=10, e =2.52 

I Phase 1 ; ġ=27r/3 


L7 - Height 7 ; H,5, 20.723 


At = 1.22m 
Ae At 
Sonic Wave | | | 


Gauges C) Q (ጋ 


ደበፎም ምች ን s 






Rubble 
ከ=3.05፲6 : 3 : Mound 
C C Q Breakwater 


10m 


Figure 3.2- Sonic Wave Gauge Locations (Fassardi 1993) 


For a small data sample of ten seconds, Figure 3.3 shows how the clean time 


series are phase shifted at each wave gauge station. Figure 3.4 shows two minutes of the 





| 
| 
| 


clean time series at Gauges 1 and 2 and the resulting incident and reflected time series 
resolved by the Fassardi (1993) algorithm. The time series at Gauge 1 is identical to 
that at Gauge 2, except for the phase shift. The incident and reflected time series 
resolved from Gauges 1 and 2 are plotted on the same scale to show the relative 
magnitude of the incident and reflected wave amplitudes. The magnitude of the wave 
amplitudes of the reflected time series are about 25% of the incident amplitudes. Figure 
3.5 shows the composite spectrum at Gauge 2 and the incident and reflected spectra 
resolved from Gauges 1 and 2 and Gauges 2 and 3,respectively. The spectra in Figure 
3.5 have been smoothed using the non-statistical technique. Box Car Averaging. The 
resolved spectra show that the energy in the reflected spectrum is about 5-10% of that in 
the incident spectrum. Breakwater reflection of about 25% observed in the time series 
results in the reflected spectral density being 5-10% of the incident spectral density, for 
each frequency component. The low frequency spike in the spectral density at 
approximately 0.04 Hz is attributed to a seiche in the wave channel caused by the start- 
up of the wave generator. The resolved spectra in Figure 3.5 are similar, indicating that 
the Fassardi (1993) algorithm produces the same spectra from different wave gauge 
pairs. At the truncation limits of 0.2140 Hz and 0.7644 Hz, the spectral density of the 
incident and reflected waves begins to diverge because the term | sin kA¢]| in the 
denominator of the equation for the amplitudes a and b in the Fassardi algorithm 
becomes small. Further examples of the resolution of incident and reflected time series 
and spectra from the output of RECANSYG for the four envelopes discussed in 


APPENDIX B are included in APPENDIX C - Time Series and Spectra. 
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Figure 3.3- Phase lag between time series at three wave gauges 
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SUMMARY AND CONCLUSIONS 


Multiple failures of existing rubble mound breakwaters designed with the Shore 
Protection Manual's (1984) methodology have indicated a need to improve the design 
procedures for rubble mound structures. Many structural and environmental 
characteristics such as armor grading, structure permeability, wave period, storm 
duration, spectral shape, and wave groupiness have secondary effects on breakwater 
stability. Irregular waves generated solely from an energy spectrum with random phases 
do not accurately model real wave trains, and the succession of wave heights or 
groupiness is an important factor in the design of rubble mound breakwaters. The 
magnitude of the energy flux given by the envelope exceedance coefficient o is one way 
of quantifying the groupiness of a wave time series. Intermediate scale wave channel 
studies at Oregon State University documented by Medina et al. (1990), Fassardi (1993), 
and Medina et al. (1994) have demonstrated that the envelope exceedance coefficient o 
may be used to correlate wave groupiness with breakwater damage. No significant 
correlation was found between the peak enhancement factor y and rubble mound 
breakwater stability. Furthermore, the constant phase shift ¢ applied to the time series 
for each wave envelope A(x,t) did not effect the armor damage. 

To calculate the envelope exceedance coefficient «o from the incident wave time 
series, either the Kimura (1985) or the Fassardi (1993) algorithm may be used to resolve 
the incident and reflected waves from closely spaced wave gauges. The Fassardi (1993) 


algorithm, used by Medina et al. (1990) in their experiments at Oregon State University 
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is different from the Kimura (1985) algorithm by the sign of the random phase angle of 
the incident time series (Eq. 3.1a & 3.8a). This sign change does not have any effect on 
the value of « or the interpretation of the results. For the frequency range given by 
Table B.1, the Fassardi (1993) algonthm effectively resolves incident and reflected wave 
time series from closely spaced wave gauges aligned in the direction of wave 
propagation. Some divergence of the incident and reflected wave spectra is observed 
near the cut off frequencies of the truncated spectra. 

Since the wave resolution algonthms presented only resolve wave components in 
the frequency range determined by the wave gauge spacing, the spectra used to simulate 
the time series for the laboratory experiments were truncated accordingly. The variance 
of the full spectra was preserved for the truncated spectra by considering the squared 
significant wave height ratio given by Eq. (B.6). For the physical simulations performed 
by Medina et al. (1990), phase spectra resulting in the highest and lowest values for the 
envelope exceedance coefficient a were chosen for each amplitude spectra. Significant 
wave heights of consecutive runs were increased in discrete increments such that the 


Stability numbers for the two armor rock sizes were equal in consecutive runs. 
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APPENDIX A: Rayleigh Distribution 


The Rayleigh probability density function (pdf) for the wave amplitudes a (= 


H/2) is given by 


A.1 
= 1 ከ ጠው a) ees 


p(a) = —exp | - 


ጪ 





2 
R 


where ap = Rayleigh parameter (Hoffman and Karst 1975) and U(a) = Heaviside step 


function. The Rayleigh cumulative distnbution function (cdf) for wave amplitudes 1s 
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where u, = average of the amplitude a; H,,, — root-mean-square wave height; and m, = 
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variance of the time series for the water surface elevation y(t). When the random 


variable a in (A.2) is scaled by the standard deviation of the time series /m, ‚then the 


following change of variables gives: 








x da 
DES pla) TE 


P(E) > eo -A |l] 
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Eqs. (4.2) may be derived from (A.1a) by solving for the generic Rayleigh parameter ۰ 
The percent or fraction of wave amplitudes greater than wave amplitude 
a, (= H,/2), may be computed from (A.1a) by 


co 


| p(a)da 


5 
| 


(A.3) 


and the natural logarithm of (A.3) gives 
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Een (A.4) 
ፎረ 


= 


where O«nxl. The average wave amplitude greater than a, may be determined from 


z, | pajas = [ ap(ayaa 


that may be integrated by parts to obtain 


H 2 
a (n) » no,y -2ln(n) * | exp - E da 


by (A.3) and (A.4). Finally, integration gives 





0 


à, = y-2In(n) + 25 oen (A.5a) 


A 


(A.5b) 


= ¥-2In(n) + 0 | 
n 


where erf( +) = the error function and erfc( -) = complementary error function 


(Kreyszig 1983). In order to relate the Rayleigh parameter œg to the average p, set n = 


| MEG. (A.5) and 


(A.6) 


and substitution into (A.1a) gives (A.2a). 





Alternatively, (A.5) may be derived directly from special functions according to 


CO co 


a, | p(a)da | ap(a)da 


A a 
a a 


T 2 3 37 
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where the generalized incomplete Gamma function T'(a,27,5,z,) (Wolfram 1991) is defined 


ER 


ከከም ۹ = az) - 1(5,2)] 


and the incomplete Gamma function [(a,z) (Wolfram 1991) is defined as 


T(a,z) - l'(a) - G (a ,z| - | t*“'exp(-t) dt 


ሯ 


H 


G(a,z) = | SS Exp (1) dt 


Substituting for [(a,Z,Z,) 


۳2 - ሃጃ 


| 


2 


t 





| E | vt exp(-t) dt 


= 
N 
ኩጋ | Lo 
ኩጋ 
دم‎ 
سس‎ 
1 
.مر‎ 
ኩጋ | دا‎ 
t2 


vt exp(-t) dt 


1 
e 





e, 


n/-In(n) « O) 


ny-In(n) + 3 eel “in(n) | 
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giving (A.5). 
To relate the Rayleigh parameter œg to Hmm, make the following change of 


variables in (À.12): 


ya? ; ዐሃ = 2805 = 2ሃሃ 03 








and 
p(a)da = p(y)dy 
giving 
d 
E 
dy 
(A.7) 
8 ሃ 
E exp | -—— | U(y) 
Der | | 


The average of a’ may now be computed from 


co 


mm | yp (x) dy 








| 
“ሙጫ 
tO 
EIS 


R 


20 
= 1 7 (ሊ.8) 


i 
ኩጋ 
R 


and the Rayleigh parameter a, may be defined by 
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(A.9) 
Substituting (A.6) and (A.9) into (A.5b) gives 
Wi 2 E era e 
H. T n 
m = ¥-In(n) لا‎ -]n (n) | (A.10b) 
B. = y-8ln(n) + ۷2۳ زرم امی‎ (A.10c) 
n 
m 


The significant wave height H, or average of the highest 1/3 wave heights H,,; may be 
computed from (A.5b) for n = 1/3 (dropping the overbar average notation) according to 








H 
B= = 2.00215 (A. 11a) 
20, ፡ 
M | (A.11b) 
5 ጋ 00215 | 


Substituting (A.11b) into (A.5b) gives 
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aed 


(A.12) 


E, 2007500215 9:2 Tn (በስው LS selte 


HIH, 


To find the mode (or most probable) value of ay, (— H4,,/2), the maximum value of 


(A.1a) occurs when 








N a. (ል.13) 
da l 

giving 

ን = Qk ; H node = 2 0 un) 
and 
E E (A.15a,b) 
K, 7 «A, 2 
H H 

E A 001) (A.15c,d) 


n Km. 


_ The mean or median of a occurs when P(a) — 0.5 in (A.1b). The natural logarithm of 


 (ለ.15) for ፻ርፌጊ) = 0.5 15 





d ned = 21n (2) (A.16) 
0 


R 
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giving 

















dn Ara 2 و‎ (A.17a,b) 
7 x BENE يو‎ 
= /8ln(2) ፡ mei = y2In(2) (A.17c,d) 


These results are summarized in Table A.1 


height 


: (2.00215)! 


EP -1/2) ۱ /In (2) /81n (2) /21n (2) 
à 3 እ ከ 8 
3 2.00215 





ean) Y dei? 
2.00215 


4.06198 | አሠ. 5.09094 27137 


5.32423 2.35924 6.67293 1.66644 





Table A.1- Rayleigh distributed wave height relationships 
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APPENDIX B: Simulation of Sea States 


In order to correctly describe wave groups in laboratory studies, the incident wave 
train must be known. An algorithm developed at Oregon State University by Fassardi 
(1993) was used to resolve the incident and reflected wave trains in a finite bandwidth 
determined by the wave gauge spacing A? (Table 3.1). One broad (y=1) and one 
narrow (y=10) truncated Goda-JONSWAP spectrum (Goda 1985) were simulated 


according to 


٢ ን ን sec ERN. 

S (f)= Cs Hg age A exPL- 1.252)" | , (B.1) 

P 
Ge 0.0624 (B.2) 

0.230+ 0.0336y- 0.185(1.9+ y)” 

۳۱/7 ۰ ۱! 
E p (B.3) 

: Ee eo 


for fun ፎ፤ ፎ fax where f, is the peak frequency, Hp is the squared significant wave height 
ratio defined in Eq. (B.6), and ፲፲, ፌጊረ 15 the significant wave height of the truncated 
- spectrum. Eq. (B.1) is different than Eq. (30) in Fassardi (1993) in that Fassardi (1993) 
| does not document the use of the constant c, defined by Eq. (B.2). The spectra were 


| truncated according to 


| f. -0.7f 


mur P 
| 


f =2.5f 
| P 


(B.4) 
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To define the truncated spectra and preserve the variance of the full spectra, Eq. 
(B.1) includes the dimensionless coefficient Hz. If the significant wave height is given in 


terms of the vanance (Table A.1) by 


H,- 4.0043/ m, (B.5) 


then the squared significant wave height ratio Hg may be defined as 


2 
H= Hau (B.6) 
2 


s,lrunc 





where H, eu is the significant wave height of the full spectra. For peak frequencies 
between 0.l and 0.5 Hz, the average dimensionless variance preserving coefficients for 
y=1l and y=10 are H2=1.0375 and H,=1.0133 respectively. 

The mean frequency of the truncated spectra defined by the first and zeroeth 


spectral moments 


e (B.7) 


was held constant at f= 0.3333 Hz such that the peak frequency of the truncated spectra 


could be computed from the frequency ratio 


(B.8) 


nn‘ 
ዘ 
tj SE 


For peak frequencies between 0.1310 0.5 Hz, the average frequency ratios were 
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f,=0.8115and fh=0.9173 for y=1 and y=10, respectively. Accordingly, the peak 
frequencies used to simulate the truncated spectra from Eq. (B.1) were f,=0.2705 Hz 
and f,=0.3058 Hz for y=1 and y=10, respectively. 

Toe E. gauge spacing defined by Figure 3.1 was 1.22m, and the water depth 
was 3.05m (Fassardi 1993). If the band pass filter limits proposed by Goda and Suzuki 
(1976) given by Eq. (3.62) are used, and the linear wave theory dispersion relationship is 


given by 


(B.9) 





the minimum and maximum resolvable frequencies ፻0፲ ልያ =1.22m are fin = 0.1299 Hz 
and f,,,—0.7589 Hz. For cases where y=1, the minimum and maximum frequencies of 
the truncated spectra are f,,,=0.1894 Hz and f,,,=0.6762 Hz. For y=10, the minimum 
and maximum frequencies of the truncated spectra are £,,,20.2140 Hz and f,,,=0.7644 
Hz. These frequency limits for the truncated spectra are reasonably close to those 
resolvable by wave gauges spatially separated by Af =1.22m. 

Different phase spectra produce time series with different wave grouping 
characteristics. Accordingly, time series were synthesized using two different phase 


spectra for each of the truncated Goda-JONSWAP spectra given by Eq. (B.1) with 
y=1,10and f=0.3333 Hz. The phase spectra were chosen from 100 random 
Deterministic Spectral Amplitude (DSA) simulations for each amplitude spectra. Values 


of aw were calculated from each DSA simulated time series according to Eq. (2.5). For 
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each y, the phase spectra that produced the wave envelope A(x,t) with the highest and 


lowest wave height variability characterized by a were chosen for the physical 


simulations. The parameters used for wave simulation are summarized in Table B.1. 












resolvable for ۱ 
ሰ መ! 
‘om Eo GOODE 


min s 


Table B.1 - Summary of simulation parameters for Af=1.22m and h=3.05m 
(frequencies in Hz) 





Low a | High a 









Realizations with the same envelope, but different phase-shifted wave profiles are 
given by (Medina et al. 1994) 
M 
n(x, t= >> a cos[kx-wt-(e+¢)] (B. 10) 
n=l 
where ¢ is a constant phase shift applied to all M wave components. Three different 
phase shifted realizations (¢=0, 27/3, 47/3) and the replicate (4=27r) for each of the 
four wave envelopes were tested to determine if realizations shifted by a constant phase 


produce different levels of damage. For each y, A(x,t), and constant phase shift $, a 


sequence of time series were tested with increasing significant wave heights. 


SS 








The breakwater model tested by Medina et al. (1990) was divided in the direction 
of wave propagation into two sides with a different armor layer rock size (viz., W, =128.5 
N and W,-99.1 N) on each side. The monochromatic wave height corresponding to the 


zero damage condition (SPM 1984) can be calculated from 


(B.11) 





where W is the median value of the mass distribution of rocks in the armor; Ky is the 
stability coefficient; O is the angle of the breakwater slope measured from horizontal, 
and p, and p, are the weight densities of the rocks and water, respectively. Solving Eq. 
(B.11) numerically for the small rock size (W,=99.1 N) yields a monochromatic wave 
height of H=0.55m for the zero damage condition. Breakwater model and armor rock 


characteristics are given in Table B.2. 


Table B.2 - Summary of breakwater model characteristics (Fassardi 1993) 





For the zero damage condition, if the representative design wave height for a 
random sea is Hno, then the design significant wave height for the small rock size is 


determined (Table A.1) by 
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ሺን > A o W.299.] N (B.12) 


۱ 





An empirically derived number which controls the stability of the armor layer (Medina et 


al. 1990) is given by 


ar SE (B.13) 


where N,is the stability number, and H, is the design significant wave height. The 
significant wave heights in the physical simulations were increased in discrete increments 
such that the stability numbers given by Eq. (B.13) for the large (W,) and small (Ws) 
rocks were equal in consecutive runs. Seven significant wave heights were tested for 
each of the sixteen realizations. The significant wave heights varied from the design 
significant wave height of the small rock size to the maximum wave height, avoiding 
breaking. From Eq. (B.13), the significant wave height of the Kk" run is given by 


۴-1 


ነከ aa 
H es aL SN ۹٩ (B.14) 
5. VW 


5 


Thus, twenty-eight realizations of four phases and seven wave heights were 
simulated for each of four target envelopes shown in Figure B.1. For each realization, 
runs of approximately 30 minutes long containing N=2'°=32768 points sampled at 


At=0.06 seconds were simulated. 
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APPENDIX C: Wave Time Series and Spectra 





Section Pace 
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MEMINI o nenne 40 
A EE 5, 06 06 59 9 65፡6 59 585549 45 
Eee eines en 50 
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RESOLVING INCIDENT AND REFLECTED TRAINS FOR RECORD: elpl!7 


Gauge No.:1 


Raw Time Series Mean Value : 0.009327 Variance : 0.040746 
Gauge No. : 2 

Raw Time Series Mean Value : 0.010709 Variance : 0.038375 
Male No. : 3 

Raw Time Series Mean Value : 0.006861 Variance : 0.036691 


RESOLVED INCIDENT AND REFLECTED WAVES 
Analysis from Gauges 2 and 3 


Incident wave Spectrum's Variance: 0.032890 
Reflected wave Spectrum’s Variance: 0.001678 


Breakwater’s reflection coefficient : 0.2259 

Incident significant wave height : 0.7261 

| Reflected significant wave height : 0.1640 

| Mean Period (TO1) from incident spectrum : 3.065 


Incident wave train, Mean and Variance from TS(2,3) 
ስክ ስ. 0.000000 TS Variance : 0.032890 


Reflected wave train, Mean and Variance from TS(2,3) 
TS Mean : 0.000000 TS Variance : 0.001678 


Analysis from Gauges 1 and 2 


Incident wave Spectrum's Variance: 0.033244 
 " ከፒ ምን wave Spectrum's Variance: 5 


Excalevater's reflection coefficient : 0.2587 

Incident significant wave height : 0.7300 

Reflected significant wave height : 0.1889 

Mean Period (TO1) from incident spectrum : 3.074 





|! Incident wave train, Mean and Variance from TS(1,2) 
Meeriean : 0.000000 TS Variance : 0.033244 


Reflected wave train, Mean and Variance from TS(1,2) 
TS Mean : 0.000000 TS Variance : 0.002225 
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RESOLVING INCIDENT AND REFLECTED TRAINS FOR RECORD: e2pll7 


Gauge No.: 1 

Raw Time Series Mean Value : 0.009919 Variance : 0.039530 
Gauge No. : 2 

Raw Time Series Mean Value : 0.003559 Variance : 0.041591 
Cee No. : 3 

Raw Time Series Mean Value : -0.010094 Variance : 0.041708 


E TED INCIDENT AND REFLECTED WAVES 
Analysis from Gauges 2 and 3 


Incident wave Spectrum's Variance: 0.031342 
Reflected wave Spectrum's Variance: 0.002484 


Breakwater’s reflection coefficient : 0.2815 
Incident significant wave height : 0.7089 

Reflected significant wave height : 0.1996 

Mean Period (TO1) from incident spectrum : 3.245 


Incident wave train, Mean and Variance from TS(2,3) 
TS Mean : 0.000000 TS Variance : 0.031342 


Reflected wave train, Mean and Variance from TS(2,3) 
TS Mean : 0.000000 TS Variance : 0.002484 


Analysis from Gauges 1 and 2 


Incident wave Spectrum's Variance: 0.033780 
Reflected wave Spectrum’s Variance: 0.002935 


Breakwater’s reflection coefficient : 0.2948 
Incident significant wave height : 0.7359 

Reflected significant wave height : 0.2169 

Mean Period (T01) from incident spectrum : 3.244 


Incident wave train, Mean and Variance from TS(1,2) 
TS Mean : 0.000000 TS Variance : 0.033780 


Reflected wave train, Mean and Variance from TS(1,2) 
TS Mean : 0.000000 TS Variance : 0.002935 
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Resolved from Gauges 2,3 (m) 


Resolved from Gauges 2,3 (m) 
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RESOLVING INCIDENT AND REFLECTED TRAINS FOR RECORD: e3p117 


Gauge No. : 1 

Raw Time Series Mean Value : 0.011743 Variance : 0.041173 
Gauee No. : 2 

Raw Time Series Mean Value : 0.002714 Variance : 0.038172 


Ganee No. : 3 
Raw Time Series Mean Value : -0.002450 Variance : 0.038333 


RESOLVED INCIDENT AND REFLECTED WAVES 
Analysis from Gauges 2 and 3 


Incident wave Spectrum’s Variance: 0.035032 
Reflected wave Spectrum’s Variance: 0.001793 


Breakwater’s reflection coefficient : 0.2262 

Incident significant wave height : 0.7494 

Reflected significant wave height : 0.1695 

Mean Period (TOl) from incident spectrum : 3.056 


Incident wave train, Mean and Variance from TS(2,3) 
TS Mean : 0.000000 TS Variance : 0.035032 


Reflected wave train, Mean and Variance from TS(2,3) 
TS Mean : 0.000000 TS Variance : 0.001793 


Analysis from Gauges 1 200 2 


Incident wave Spectrum's Variance: 0.034591 
Reflected wave Spectrum’s Variance: 0.002552 


a Water s reflection coefficient : 0,2716 

Incident significant wave height : 0.7447 

Reflected significant wave height : 0.2023 

Mean Period (T01) from incident spectrum : 3.078 


Incident wave train, Mean and Variance from TS(1,2) 
TS Mean : 0.000000 TS Variance : 0.034591 


Reflected wave train, Mean and Variance from TS(1,2) 
TS Mean : 0.000000 TS Variance : 0.002552 
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Resolved from Gauges 2.3 (m) 


Resolved from Gauges 2,3 (m) 
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BESSIEVING INCIDENT AND REFLECTED TRAINS FOR RECORD: 7 


Gauge No. : } 
Raw Time Series Mean Value : 0.015713  Vanance : 0.043060 


Gauge No.:2 
Eme Senes Mean Value : 0.02582] Variance : 0.041359 


Gauce No. : 3 
Raw Time Series Mean Value : 0.029771 Variance : 07 


E VED INCIDENT AND REFLECTED WAVES 
Analysis from Gauges 2 and 3 


Incident wave Spectrum's Variance: 0.036691 
Reflected wave Spectrum's Variance: 0.002226 


Ereakwater's reflection coefricient : 0.2463 

Incident significant wave height : 0.7670 

Reflected significant wave height : 0.1889 

Mean Period (TOl) from incident spectrum : 3.206 


Incident wave train, Mean and Variance from TS(2,3) 
TS Mean : 0.000000 TS Variance : 0.036691 


Reflected wave train, Mean and Variance from TS(2,3) 
۲۱۱۱ 62۲۱ : 0.000000 TS Variance : 0.002226 


Analysis from Gauges 1 and 2 


Incident wave Spectrum's Variance: 0.036596 
Eeilecied wave Spectrum's Variance: 0.002472 


Breakwater’s reflection coefficient : 0.2599 
Incident significant wave height : 0.7660 

Reflected significant wave height : 0.1991 

Mean Period (101) from incident spectrum : 3.212 


Incident wave train, Mean and Variance from TS(1,2) 
TS Mean : 0.000000 TS Variance : 0.036595 


Reflected wave train, Mean and Variance from TS(1,2) 
TS Mean : 0.000000 TS Variance : 0.002472 
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Resolved from Gauges 1,2 (m2/H2) 


Resolved from Gauoes 2,3 (m 2/2) 


Spectral Density 0ቹ፳1 (8111) 


Spectral Density a Gauge 2 (m21) 
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APPENDIX E: Notation 


wave envelope function 
FFT coefficient at i" wave gauge 
the median of a 


most probable value of a 


the average wave amplitude greater than a, 


root mean square wave amplitude 

incident wave amplitude 

FFT coefficient at i" wave gauge 

reflected wave amplitude 

Goda-JONSWAP spectrum parameter defined in Eq. (B.2) 
wave amplitude interval 

wave height interval 

interval of & defined in Eq. (A.2e) 


expected value of a’ 


EROS UNCON Of (°) 


complementary error function of (*) 


wave frequency 


mean frequency of spectrum 


— maximum (minimum) cut-off frequency 
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H(t) 
111) 
H(nAt) 


med 


Ha. 


= peak frequency of spectrum 


= frequency ratio given by Eq. (B.8) 


eamma function from APP. A 

gravitational constant 

wave height 

wave height function at a fixed location 

wave height function 

discrete wave height function at x=0 

the median of H 

the most probable value of H 

Squared significant wave height ratio given by Eq. (B.5) 
root mean square wave height 

significant wave height = Hi» 

significant wave height of full spectrum 

significant wave height of the k” run 

significant wave height of truncated spectrum 

average of the wave heights 

average of the highest one-third of the wave heights 
average of the highest one-tenth of the wave heights 
average of the highest one-hundredth of the wave heights 


water depth 
ገ 
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Kp 
k (=2r/L) 
E 


nin) 
mo 


m; 


EC 
D: 
5,() 


armor layer stability coefficient 

wave number 

wave length 

maximum (minimum) wave length 

variance of the time series, zero-th spectral moment 
first spectral moment 

percent or fraction of wave amplitudes greater than a, 
number of data points 


stability number defined ın Eq. (B.12) 

Rayleigh cumulative distribution function of (°) 
Rayleigh probability density function of ( ፥ ን 
Goda-JONSWAP spectrum defined by Eq. (B.1) 
wave period 

time 

Heaviside step function of (-) 

the median value of the mass distribution of rocks in the armor 
large armor layer rock size 

small armor layer rock size 

complex valued FFT coefficient 

position of i^ wave gauge from wave maker 
change of variable from Eq. (A.7) 


envelope exceedance coefficient 
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Pa 
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Py 
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Yim 


w(=2r/T) 


magnitude of the variation of wave energy above Hun 
Rayleigh parameter 

reflected wave phase 

incomplete gamma function from APP. A 

generalized incomplete gamma function from APP. A 
peak enhancement factor 

normalized measure of the variation of wave height 
wave gauge spacing 

Sampling time interval 

incident wave phase 

reflected time series at i" wave gauge 

sea surface elevation 

incident time series at i^ wave gauge 

angle of the breakwater slope measured from horizontal 
average of the iude a 

change of variable from Eq. (A.2e) 

weight density of rock 

weight density of water 

Goda-JONSWAP spectrum parameter defined in Eq. (B.3) 
constant phase shift applied to each realization 
incident (reflected) spatial wave phase 


circular wave frequency 
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